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PREFACE 


The  Electromagnetic  Compatibility  Analysis  Center  (ECAC)  is  a Department  of 
Defense  facility,  established  to  provide  advice  and  assistance  on  electromagnetic 
compatibility  matters  to  the  Secretary  of  Defense,  the  Joint  Chiefs  of  Staff,  the  military 
departments  and  other  DoD  components.  The  center,  located  at  North  Severn,  Annapolis, 
Maryland  21402,  is  under  the  policy  control  of  the  Assistant  Secretary  of  Defense  for 
Communication,  Command,  Control,  and  Intelligence  and  the  Chairman,  Joint  Chiefs  of 
Staff,  or  their  designees,  who  jointly  provide  policy  guidance,  assign  projects,  and  establish 
priorities.  ECAC  functions  under  the  executive  direction  of  the  Secretary  of  the  Air  Force 
and  the  management  and  technical  direction  of  the  Center  are  provided  by  military  and  civil 
service  personnel.  The  technical  operations  function  is  provided  through  an  Air 
Force-sponsored  contract  with  the  I IT  Research  Institute  (IITRI). 

This  report  was  prepared  as  part  of  AF  Project  649E  under  Contract 
F-I9628-78-C-0006  by  the  staff  of  the  NT  Research  Institute  at  the  Department  of  Defense 
Electromagnetic  Compatibility  Analysis  Center. 

To  the  extent  possible,  all  abbreviations  and  symbols  used  in  this  report  are  taken  from 
American  Standard  Y10.19  (1967)  "Units  Used  in  Electrical  Science  and  Electrical 
Engineering"  issued  by  the  USA  Standards  Institute. 

Users  of  this  report  are  invited  to  submit  comments  that  would  be  useful  in  revising  or 
adding  to  this  material  to  the  Director,  ECAC,  North  Severn,  Annapolis,  Maryland  21402, 
Attention  XM. 
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EXECUTIVE  SUMMARY 


Computer  models,  developed  using  the  Taylor- Fourier  series 
method  to  analyze  double-balanced-diode  and  single-diode  mixers, 
are  described.  Use  of  this  method  results  in  mixer  models  that 
are  considerably  different  in  form  than  those  obtained  by  past 
methods;  the  predicted  response  levels  are  in  close  agreement 
with  measured  data.  These  models  were  developed  by  the  Electro- 
magnetic Compatibility  Analysis  Center  (ECAC)  to  provide  improved 
procedures  for  estimating  desired  and  undesired  receiver  responses 
due  to  the  nonlinearity  of  the  mixer.  These  undesired  responses 
may  include  cross  modulation,  intermodulation,  spurious  responses 
and  desensitization. 

The  models  are  used  to  compute  nonlinear  response  amplitudes 
based  on  the  nominal  values  of  resistors,  inductors,  and  capaci- 
tors in  the  circuit,  the  local  oscillator  voltages,  and  the 
direct  current  (dc)  characteristic  of  the  diode (s).  Measured 
spurious  responses  were  compared  with  spurious  responses  calcu- 
lated by  the  models  and  were  in  close  agreement. 

Procedures  for  developing  additional  mixer  models  are  also 
described. 

The  methods  described  herein  and  the  measured  intermodulation 
data  normally  acquired  for  equipment  in  obtaining  a spectrum  sig- 
nature can  be  used  to  make  reasonable  estimates  of  levels  of  many 
other  intermodulation  responses.  These  other  responses  can  be 
significant  in  an  electromagnetic  compatibility  analysis,  but  are 
not  measured. 
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SECTION  1 
INTRODUCTION 


BACKGROUND 


Frequently,  EMC  analyses  require  estimation  of  the  effects  of  non- 
linear interactions.  Equipment  proximity  and  the  resulting  high  inter- 
fering power  levels  on  platforms  such  as  airplanes,  ships,  and  satellites, 
and  in  tactical  environments,  may  produce  the  following  nonlinear  inter- 
actions : 

1.  Cross  modulation 

2.  Receiver  intermodulation  (RIM) 

3.  Spurious  responses  (SR) 

4.  Desensitization 

5.  Transmitter  intermodulation  (TIM). 

Estimates  are  sometimes  based  on  measured  data  taken  on  a small 
sample  of  equipments  of  the  same  or  similar  design.  However,  because 
of  the  limited  scope  and  increasing  costs  of  most  measured  data,  com- 
puterized equipment  models  are  becoming  more  prevalent. 

In  the  past,  a major  stumbling  block  in  the  development  of  accurate 
nonlinear-interaction  receiver  models  has  been  the  modeling  of  the  mixer,1 
in  many  cases  because  of  a basic  flaw  in  the  modeling  method.  A commonly 
used  method  can  be  described  as  follows.  The  current-voltage  transfer 


Donaldson,  Jr.,  E.  E.  and  Moss,  R.  W. , Study  of  Receiver  Mixer  Charac- 
teristics , Final  Report,  Technical  Report  ECOM-01426-F,  USAEC,  Fort 
Monmouth,  NJ,  September  1966. 
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characteristic  is  represented  by  a power  series;  the  input  to  the 
mixer  is  the  sum  of  the  radio-frequency  (RF)  and  local  oscillator 
(LO)  signals.  For  practical  mixer  design,  a flaw  exists  because 
the  LO  level  is  so  high  that  either  the  series  does  not  converge 
or  the  required  number  of  terms  becomes  prohibitive  (see  Section 
2 and  APPENDIX  G) . 

OBJECTIVE 

* The  objective  of  this  analysis  was  to  present  improved  proce- 

dures for  estimating  receiver  responses  resulting  from  mixer  non- 

► 

linearities . 

APPROACH 

A literature  search  was  conducted  for  available  mixer  analysis 
methods.  Most  of  the  methods  found  were  either  not  applicable  to 
common  mixer  circuits  or  were  not  validated.  One  of  the  nonvali- 
dated  methods2’3  that  appeared  promising  was  selected  for  validation. 

This  method,  which  uses  the  Taylor-Fourier  Series  mixer  model, 
is  described  in  Section  3.  It  is  used  to  calculate  the  magnitudes 
of  response  coefficients  of  a mixer.  These  response  coefficients 
may  be  used  to  estimate  the  levels  of  the  mixer's  desired  and  undesired 
responses  (i.e.,  spurious  responses,  intermodulation  responses). 

A 5-step  method  was  used  for  developing  the  Taylor-Fourier  series 
model  for  mixers  and  is  presented  in  Section  3.  Mixers  are  modeled 
as  memoryless,  voltage-dependent,  current  sources.  The  steps  were 
as  follows. 


2Ehrman,  L.  and  Graham,  J.  eds..  Nonlinear  System  Modeling  and  Analysis 
with  Applications  to  Communications  Receivers , RADC-TR-73-178 , USAF 
Rome  Air  Development  Center,  Rome,  NY,  June  1973. 

3Terman,  F.  E.,  Electronic  and  Radio  Engineering,  McGraw-Hill,  4th  Ed., 
New  York,  NY,  1955.  Pp  573-575. 
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1.  Model  the  current  versus  voltage  characteristic  of 
the  nonlinear  device. 

2.  Write  the  loop  or  node  equations  for  the  circuit. 

3.  Derive  the  expressions  for  the  derivatives  of  the 
load  current  with  respect  to  the  RF  source  voltage.  (A  numerical 
procedure  could  be  employed.) 

4.  Derive  the  expressions  for  the  Taylor  series  coefficients 
as  time  variables. 

5.  Calculate  the  Taylor-Fourier  series  coefficients  from 
a Fourier  series  expansion  of  each  Taylor  serie-s  coefficient. 

Computer  routines  based  on  the  5 -step  method  were  then  written  and 
programmed.  The  computer  routines  were  exercised  to  determine  the 
theoretical  response  coefficient  levels  for  a single-diode  mixer  (SDM) 
and  a Hewlett  Packard  double-balanced-diode  mixer  (DBDM)  circuit. 

To  help  validate  the  Taylor-Fourier  Series  mixer  models,  measurements 
were  made  at  ECAC  of  the  mixer  spurious  responses  resulting  from  the  SDM 
and  the  DBDM  circuits.  Comparisons  of  the  model  predictions  and  the 
measured  data  were  performed  and  are  presented  in  Section  4 and  APPENDIX 
C. 


Procedures  for  calculating  the  receiver  spurious  response  rejection 
on  the  basis  of  the  mixer  spurious-response  coefficients  are  described 
in  APPENDIX  B. 


A procedure  for  estimating  certain  high-order  RIM  coefficients  on 
the  basis  of  measured  values  of  low -order  RIM  coefficients  is  described 
in  APPENDIX  A.  The  rationale  for  this  extrapolation  is  based  on  the 
results  of  mixer  performance  predictions  utilizing  the  Fourier-Series 
model. 
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Finally,  a comparison  of  the  modeling  differences  between  the 
constant  coefficient  power  series  and  the  Taylor-Fourier  Series  is 
presented  in  Section  5. 
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SECTION  2 

THE  POWER  SERIES  MIXER  MODEL 

The  power  series  analysis  technique  has  been  used  to  describe  the 
spurious  and  intermodulation  response  characteristics  of  a mixer.  The 
results  are  considerably  different  (Section  5)  from  those  obtained  by 
the  Taylor-Fourier  series  analysis  technique,  as  described  in  Section  3. 
As  indicated  below,  the  power  series  is  not  adequate  when  applied  to  a 
practical  mixer  design.  The  reason  for  this  is  that  the  power  series 
coefficients  are  not  constant  with  time  due  to  the  fact  that  the  LO 
level  is  so  large  in  practice  that  the  "operating  point  level"  (i.e., 
where  the  derivatives  are  evaluated)  changes  significantly  over  the 
time  of  one  period  of  the  LO  waveform.  Yet  in  a power  series  analysis, 
the  coefficients  are  used  as  time-constant  values.  The  Taylor-Fourier 
series  takes  the  time  variation  into  account  by  expanding  each  Taylor 
(or  power)  series  coefficient  in  a Fourier  series. 

POWER  SERIES  ANALYSIS  FOR  A DIODE  MIXER 

Diode  mixer  circuits  (Figure  1)  have  been  analyzed  in  other  reports 
and  publications  (see  References  1,4,  and  5).  Those  results  are  useful 
for  the  purposes  of  this  discussion.  The  relationship  between  the  vol- 
tage (v)  and  the  current  (i)  is  given  using:5 6 

* ■ i0  [•  “<V  ' iRT)  - l]  (1) 


4Steiner,  J.  W.,  "An  Analysis  of  Radio  Frequency  in  Interference  Due  to 
Mixer  Intermodulation  Products,"  IEEE  Transactions  on  EMC,  Vol.  EMC-6, 
No.  1,  January  1964. 

5Herishen,  J.  T.,  "Diode  Mixer  Coefficients  for  Spurious  Response  Predic- 

tion," IEEE  Transactions  on  EMC,  Vol.  EMC-10,  No.  4,  December  1968. 

6Ebers  J.  J.  and  Moll,  J.  L.,  "Large-Signal  Behavior  of  Junction  Transis- 
tors," Proceedings  of  the  IRE,  December  1954. 
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where 


iQ  = reverse  leakage  current,  a constant,  in  amperes 

R^,  = the  total  series  circuit  resistance,  including  the 
diode  bulk  resistance,  R„,  in  ohms 

D 

a = a constant  whose  value  generally  lies  between  10  and 
40  volts  1 


e = base  of  the  natural  logarithm  = 2.718. 


Figure  1.  Diode  mixer  circuit. 


The  power  series  representation  of  Equation  1 is: 


i 

■ L k»c*-vb)n. 

n 

where 

It 

1 3ni 

n 

n!'  3vn 

v=Vb 

vb 

= the  bias  voltage 

I 
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In  Figure  2,  the  magnitude  of  the  power  series  coefficients  rela- 
tive to  the  second-order  coefficient  (in  dB)  are  shown  for  a diode 
mixer  (Reference  5).  The  abscissa  value,  i,  is  the  current  that  flows 
for  a given  voltage,  v. 

Based  on  a fixed  value  of  mixer  forward  bias  voltage,  a value  of 
bias  current  is  obtained  using  the  graph  in  Figure  3.  (Figure  3 was 
obtained  using  Equation  1.)  The  normalized  power  series  coefficients 
are  obtained  using  the  graphs  in  Figure  2.  The  coefficients  are  then 
substituted  into  Equation  2 from  which  the  nonlinear  interactions  are 
calculated. 


CUMMfMT.  . (»A> 
(•I  «'  S.4  AMO  9 


CUftftCMT.  . (pA) 
(Ih'lAMT 


Figure  2.  For  diode  mixer  power  series  coefficients 

normalized  to  second-order  coefficient,  in  dB. 
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AN  ILLUSTRATIVE  EXAMPLE  OF  POWER  SERIES  LIMITATION 


The  following  example  illustrates  the  inapplicability  of  the  power 
series  for  diode  mixer  analysis.  Assume  that  the  LO  voltage  is  a sinu- 
soid of  peak  amplitude  equal  to  0.63  volts.  (This  would  correspona  to 
an  available  power  of  0 dBm  for  a 5C-ohm  signal  generator.  Virtually 
all  mixers  employ  this  level  or  greater  to  provide  low  conversion  loss 
to  the  desired  signal.)  From  Figure  3 it  is  seen  that  the  current  flow- 
ing at  the  positive  peak  of  the  LO  waveform  is  approximately  2.5  mA 
(2,500  pA),  five  times  the  maximum  abscissa  value  on  Figure  2.  In  fact, 
the  abscissa  value  will  vary  from  -i  to  2,500  pA  over  each  period  of 
the  LO  waveform.  In  general,  V^  (the  true  dc  bias)  is  small  compared 
to  the  peak-to-peak  variation  in  the  LO  voltage.  Thus,  intuitively  it 
is  expected  that  one  particular  value  of  20  log  k2/kn  obtained  from 
Figure  2 will  not  suffice.  The  coefficients  will  periodically  take  on 
all  values  shown  on,  and  many  more  off,  the  scale  of  the  figure.  It  is 
difficult  to  define  a precise  voltage  limit  for  the  usefulness  of  the 
power  series.  The  number  of  terms  needed  for  adequate  convergence 
increases  with  increasing  ac  voltage  levels  and  is  dependent  on  the  bias 
current.  An  optimistic  limit  is  stated  in  the  following  inequality: 


LO 


(3) 


where 


iL0  = the  peak  of  the  LO  current  waveform. 

In  the  previous  example,  a = 20  volts'1,  R^,  = 100  ohms,  and  iLQ  < 0.5  mA 
will  satisfy  the  inequality.  The  0.5  mA  corresponds  to  an  LO  power  of 
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-14  dBm.  An  LO  power  of  0 dBm  corresponds  to  2.5  mA.  For  a = 40  volts' 
and  R^,  = 500  ohms,  the  LO  power  would  be  limited  to  -34  dBm  by  the  pre- 
ceeding  inequality.  Thus,  the  power  series  is  not  useful  for  typical  LO 
power  levels. 

It  is  shown  in  APPENDIX  G that  the  power  series  for  vacuum  tubes  will 
not  converge  in  some  practical  situations. 


10 
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SECTION  3 

TAYLOR-FOURIER  SERIES  MIXER  MODEL 


BACKGROUND 


A literature  review  of  mixer  analysis  methods  was  conducted.  A 
promising  method  was  noted  in  References  2 and  3,  but  quantitative 
results  were  not  presented. 

The  primary  characteristic  of  the  Taylor- Fourier  method  is  the 
use  of  the  Taylor-Fourier  series.  With  this  method,  it  is  assumed 
that  the  LO  level  is  much  greater  than  the  RF  input  levels.  (See 
APPENDIX  C for  additional  discussion  of  this  point.)  This  method 
was  tested  as  part  of  this  task  to  see  if  it  could  be  successfully 
applied  to  practical  mixer  analyses  and  could  provide  an  accurate 
physical  model  of  a mixer. 

The  Taylor-Fourier  method  was  applied  to  two  types  of  diode 
mixers:  single-diode  mixers  (SDM)  and  double-balanced-diode  mixers 
(DBDM) . Computer  models  for  predicting  spurious  responses  based 
on  circuit  diagrams  were  prepared.  For  both  types  of  mixers,  the 
predictions  were  based  only  on  the  circuit  passive  component  values, 
the  LO  voltages,  a reference  output  voltage,  and  the  diode  charac- 
teristics. Mixer  spurious -response  measurements  were  made  at  ECAC 
and  compared  with  predicted  responses  (Section  4).  The  procedures 
used  to  estimate  the  diode  parameters  on  the  basis  of  static  or  dc 
measurements  are  described  in  APPENDIX  D. 

THE  TAYLOR-FOURIER  METHOD 


The  procedures  employed  to  model  the  SDM  and  DBDM  apply  in 
general  to  any  mixer  represented  as  a resistive,  voltage-dependent, 
current  source.  The  output  current  is  related  to  the  input  RF  vol- 
tage, vRp, using: 
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1 ■ Jj  (”RF>Q  aP,Q  c°s  <P2*fl.0t)  <4> 


where 


p.Q 


integers  representing  the  harmonic  number  of  the  LO 
and  RF  inputs,  respectively,  contributing  to  the 
response  being  analyzed 

the  Taylor-Fourier  series  response  coefficients 

the  local  oscillator  frequency 

time. 


USING  THE  RESPONSE  COEFFICIENTS, 


Once  the  response  coefficients  are  known,  the  mixer  output  levels 

may  be  calculated  from  Equation  4 for  any  desired  or  undesired  input 

voltage.  However,  Equation  4 can  be  put  in  a more  convenient  form  for 

calculation  purposes  if  the  form  of  v is  given.  For  example,  suppose 

Kr 

Vjjp  consists  of  a single  tone  at  a frequency  that  will  cause  a mixer 
response.  Let: 

VRF  = VRF  C0S  27'fRFt  (5) 


where 


Vrf  = peak  voltage 
fRF  = radio  frequency. 

A spurious  response  is  said  to  occur  when: 
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RF 


Pf  + f 
^rLO 

Q 


IF 


(6) 


where 


IF  = receiver  intermediate  frequency 
LO  = receiver  local  oscillator  frequency. 


Note,  the  desired  response  is  among  the  set  of  responses  described 
by  Equation  6,  but  cannot  be  termed  spurious.  If  the  IF  is  the  differ- 
ence between  the  RF  and  LO  frequency  and  if  the  LO  is  below  the  receiver's 
tuned  frequency,  then  Equation  6 reduces  to: 

fDC  = ftri  + fT_  receiver  tuned  frequency  (7) 

Kr  LU  I r , 

where 

P = Q = 1. 


In  addition,  the  term  of  interest  in  Equation  4 that  yields  a com- 
ponent at  fjp  when  the  input  signal  frequency  satisfies  Equation  7 is: 

iD  4 (VRF  C0S  2*fRFtJ  ai,l  C0S  2”fLOt 

- al,lV2  c0»t2*(*RF-  V11*  •••  (8) 

where 

i = the  desired  output  of  the  mixer  and  additional 
D 

outputs . 
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Thus,  a.  /2  is  the  conversion  transconductance  of  the  mixer. 

(For  this  case,  the  P=l,  Q=3  term  of  Equation  4 will  also  yield  a com- 
ponent at  fTC.  In  fact,  it  is  this  term  which  describes  the  dynamic 
1 r 

range  of  the  mixer.) 


The  more  general  case,  where  fRp  is  not  restricted  by  Equation  7, 
may  take  on  any  value  which  satisfies  Equation  6 for  any  integers  P and 
Q,  and  either  sign.  For  example,  setting  P=2  and  Q=3  and  sign  = +,  an 
input  of  frequency 


2fL0+  fIF 


RF 


(9) 


will  cause  a response  at  the  mixer  output  at  the  IF  stage  tuned  fre- 
quency, fjp.  The  amplitude  of  this  response  may  be  calculated  by  con- 
sidering the  appropriate,  P=2,  Q=3  term  of  Equation  4: 

iSPUR  4 (VRF  C0S  2"fRFt)3  “2.3  C0S  4,V  <l0a> 

= M'’V3a2,3  cos  [2'<2fLO-  ^rfW*---  U°>>> 

Using  Equation  9: 

1SPUR  = (Ji)3VRF3  a2  3 cos  27rfIpt  ....  (10c) 

ThUS  a2.3  VRF/8  is  the  conversion  transconductance  of  the  mixer 
for  a P=2,  Q=3  spurious  response.  A more  general  expression  for  Equa- 
tion 10c  is: 
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i - V ^ 

SPUR  " ^ RF 


aP,Q  C0S 


2-nf  T_t 
IF 


P >0 


(11a) 


and 


ispuR  = VRFQaP,Qcos  P=0  (llb) 

Equation  21  of  APPENDIX  C presents  this  relationship  in  terms 
of  power  output  (referred  to  the  input)  versus  power  available  at 
the  input.  Similar  expressions  may  be  employed  to  cover  2-signal 
intermodulation  and  other  nonlinear  effects. 

CALCULATING  THE  RESPONSE  COEFFICIENTS 


In  order  to  calculate  the  a p ^ values  the  following  five  steps  were 
used  and  are  suggested  for  use  in  developing  additional  models. 

1.  Model  the  current  versus  voltage  characteristics  of  the 
nonlinear  device.  Figure  4 i an  equivalent  circuit  diagram  of  the  Ebers - 
Moll  diode  model  (see  Reference  6).  The  nonlinearity  is  represented  by 
the  relationship  between  v^  and  i^  using: 


II 

Q 

r av' 
D 

-0 

(12) 

where 

XD 

current 

through  the  ideal 

junction,  in  amperes 

xo 

reverse 

leakage  current, 

in  amperes 

VD  = 

voltage 

across  the  ideal 

diode  junction,  in  volts 
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a = diode  parameter  = q/kT,  volt” 

q = electric  charge  of  an  electron 

k = Boltzman's  constant 

T = temperature  in  °K. 

The  value  of  the  shunt  resistance  is  generally  very  large,  and 

therefore  ignored  in  the  mixer  models.  A computer  program  was  prepared 
to  calculate  the  values  a and  R^  that  fit  a set  of  measured  dc  character- 
istic data,  as  described  in  APPENDIX  D. 


DIODE 


cl 

An  ideal  diode  is  governed  by  Equation  12. 
Figure  4.  Ebers-Moll  diode  model. 
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where 

R + R,  + R.  = total  resistance 
s b L 

RF  source  resistance 
bulk  resistance 
load  resistance 

Step  2 was  also  accomplished  for  the  DBDM.  Figure  6 is  the 
circuit  diagram  for  the  DBDM  mixer.  For  this  circuit,  four  loop 
equations  are  necessary;  ideal  transformers  are  assumed.  These 
equations  are  given  in  matrix  notation  by: 


Rt 

R 

s 

*b  " 


NOTE : 


’ nO 


T 


Eo 


D1  = diode  #1 
D2  = diode  #2 
D3  = diode  #3 
D4  = diode  #4 


Figure  6.  Modeled  DBDM  circuit. 
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where 


I in/h  + A + k3E0-kies 


- In  *2  + 1 


1 In  ( x3  + 1 


°3 


1 In  / 4 + 1 


°4  \104 


- k3Vk2es 


♦ k4E0+kl6s 


- k.En+k-e 
4 0 2 s 


• H H 


k,  = n /n 
1 ae  s 


k^  = indicated  turns  ratio,  secondary  to  primary-  k2  = n^/n^ 


Eq  = instantaneous  LO  voltage 


k3  = nbc/n0 


instantaneous  RF  voltage 


k4  = ncf/n0 


Rq  = LO  source  resistance 

Rg  = RF  source  resistance 

i^  = reverse  leakage  current  for  diode 

cr  = diode  parameter  for  diode 

**bj  = bulk  resistance  for  diode. 


The  members  of  [i]  are  i^  where  i^  is  the  forward  current  through  j 
diode.  The  members  of  the  [A]  are  given  by: 


Al,l  • -k2  Rs  - k2.  R0  - *L  - R,,, 


A.  , ■ -k.k.R  + k,  Rn  R. 

1,2  12s  30  L 
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1,3 

= 

2 

+k7R  - 
1 s 

k3k4R0 

1,4 

= 

+klk2Rs 

+ k3k4Ri 

'2,2 

= 

-kjR.  - 
2 s 

k3R0  * 1 

‘2,3 

« 

+k,k_R 

12s 

+ k3k4Ri 

‘2,4 

» 

*k2Rs  - 

k3k4R0 

‘3,3 

= 

-klRs  - 

k4R0  * 1 

‘3,4 

= 

’Ws 

* k4R0  ' 

‘4,4 

= 

-k~R  - 
2 s 

k4R0  - 1 

j , k 

s 

Ak,j  for  ^k’  ^ 

) - rl 
lL  - Rb2 
) - rl 
’ rl 

lL  ' ^3 

■ rl 

' Rb4 


3.  Derive  expressions  for  the  derivatives  of  the  load  resis- 
tor current  with  respect  to  the  RF  source  voltage.  This  was  done  for 
as  many  as  seven  derivatives  in  both  models.  The  derivatives  for  the 
SDM  are  contained  in  APPENDIX  C;  the  first  derivative  for  the  DBDM  is 
derived  using: 


where 


"Wl 

-k2 a2^2 

+k,a_J_ 
1 j 3 

+k2a4J4 


3i 
3 e 


(16) 


J. 

J 

B. 

J,J 

B.  . 
J.k 


(ij  + Vj3 

A.  .a.J.-l 
J.J  J J 

A . . a . J . 

J,k  j j 


j*k 
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The  expressions  for  the  derivative  have  the  general  form: 

[c«]  = [b 

The  expressions  for  thejc^jare  in  APPENDIX  F. 

The  load  resistor  current  is  given  using: 
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2n 

aP,Q  = f aQ(t)  C0S  (PaiL0t}  d(wL0t)  P>0  (22) 

0 

where 

WL0  = 2irf  LO 

These  coefficients  are  obtained  by  sampling  the  LO  voltage  time 
waveform.  For  each  sample,  the  current  and  a^  values  are  calculated; 
numerical  integration  is  employed.  To  use  Equation  22,  the  LO  wave- 
form must  be  an  even  function. 

FREQUENCY-DEPENDENT  RESPONSE  COEFFICIENTS 

In  some  cases,  those  frequency-dependent  circuit  parameters  which 
cannot  easily  be  isolated  may  become  significant  for  analysis  purposes. 
Some  of  these  circuit  parameters  can  be  considered  in  the  same  manner 
as  a bypass  capacitor  (a  short  circuit  at  high  frequencies  and  an  open 
circuit  at  lower  frequencies).  Where  this  is  not  possible,  the  Taylor 
series  expansion  in  the  Taylor-Fourier  series  method  must  be  replaced 
by  a Volterra  series  expansion.  The  additional  complexities  involved 
from  this  requirement  will  not  be  discussed  here. 
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SECTION  4 

RESULTS  OF  TAYLOR-FOURIER  MODELING 

Agreement  is  generally  good  between  measured  mixer  spurious  response 
levels  and  those  predicted  by  the  computerized  Taylor-Fourier  series 
models.  The  results  were  better  for  the  SDM  than  for  the  DBDM  apparently 
because  of  the  sensitivity  of  certain  responses  of  the  DBDM  to  the  dc 
diode  and  other  circuit  parameters  (especially  their  degree  of  balance). 
Other  DBDM  responses  (described  in  this  section)  are  generally  insensi- 
tive. 


The  Taylor-Fourier  models  are  described  in  Section  3 and  APPENDIX  C. 
Results  of  comparisons  of  model  predictions  with  measured  spurious 
response  (SR)  data  are  given  below.  The  predictions  were  based  only  on 
the  nominal  circuit  passive  component  values,  LO  voltages, and  dc  diode 
characteristics  (dc  current  versus  dc  voltage). 

SINGLE-DIODE  MIXERS  (SDM) 

The  measured  SR  data  for  a SDM  are  shown  in  TABLE  1(a)  of  APPENDIX 
C;7  the  corresponding  predictions  are  shown  in  TABLE  1(b).  Forty-two 
of  the  predicted  values  (over  half)  differ  by  less  than  1 dB  from  the 
measured  values.  The  average  difference  is  1 dB  and  the  standard  devia- 
tion of  the  differences  is  3.1  dB.  Some  of  the  data  are  shown  on 
Figure  7. 

A significant  conclusion  from  a parameter  sensitivity  study  was 
that  good  prediction  accuracies  are  obtainable  with  uncertainties  of 
up  to  25%  in  the  diode  parameter  a and  a factor  of  10  in  the  leakage 

7ECAC  lab  memo  *7769  May  1977. 
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q = harmonic  number  of  input  signal  frequency 

p = harmonic  number  of  LO  frequency 
(see  .Equation  C-2  of  APPENDIX  C) 


Figure  7.  Measured  and  predicted  spurious  response  values 
for  SDM. 
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current.  This  information  is  of  value  to  a desk-bound  analyst.  The 
sensitivity  analysis  results  are  presented  in  TABLES  2(a)  and  2(b)  of 
APPENDIX  C.  Note  that  the  standard  error  is  minimum  when  the  values 
of  ig  and  a,  determined  by  examining  the  dc  diode  characteristic,  are 
employed.  Similar  results  were  obtained  when  the  LO  voltage  was 
increased  from  0.47  to  1.66  volts,  as  shown  in  TABLES  4(a)  and  4(b)  of 
APPENDIX  C. 

DOUBLE- BALANCED-DIODE  MIXERS  (DBDM) 

Measured  SR  data  for  Hewlett-Packard  HP10534a  and  HP10514a  mixers 
were  obtained  at  ECAC8  as  noted  in  TABLES  1 and  2.  The  HP10534a  was 
found  to  have  a defective  diode  prior  to  testing.  Its  diodes  were 
replaced  by  four  new  1N82  diodes  after  the  dc  characteristics  of  these 
diodes  were  measured.  The  circuit  diagram  for  the  HP10534a  is  given  in 
Figure  6;  a large  number  of  DBDM's  use  this  circuit. 

DBDM  Predicted  Responses 

The  predicted  responses  for  the  HP10534a  circuit  are  given  in 
TABLE  3.  For  several  responses,  the  predictions  varied  greatly  from 
the  measured  data.  A parameter  sensitivity  study  indicated  that  these 
responses  were  among  several  that  were  sensitive  to  the  diode  para- 
meters and/or  to  the  balance  among  the  diodes. 

DBDM  Parameter  Sensitivity  Study 

The  DBDM  model  was  employed  to  perform  a study  in  which  the  diode 
parameters  --  a,  ig,  and  R^  --  were  varied  about  the  center  values  selec- 
ted by  the  procedure  of  APPENDIX  D.  The  variations  were  as  follows: 


8ECAC  lab  memo  #77-73. 
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Most  responses  remained  relatively  unchanged  for  the  input  values. 

A few  responses  were  significantly  weaker  for  a given  set  of  diode  para- 
meters. These  were  among  the  responses  that  disappear  for  a perfectly 
balanced  DBDM.  The  results  also  indicated  that  the  best  overall  predic- 
tions were  obtained  by  increasing  the  a value  of  diode  # 2 (D2,  Figure  6) 

This  change  made  diode  #2  more  nearly  equal  to  the  other  diodes 
in  parameter  values.  It  was  noted  that  the  measured  dc  data  for  this 
diode  resulted  in  the  worst  "fit"in  terms  of  rms  deviation  to  the  Ebers- 
Moll  diode  model.  The  predicted  SR  values,  witn  the  increase  of  by 
5 volts  , are  given  in  TABLE  4.  Apparently,  the  representation  of  this 
diode  was  lacking. 

A Perfectly  Balanced  DBDM 

Three-fourths  of  the  responses  (i.e.  if  either  P or  Q is  even)  com- 
pletely disappear  under  the  unlikely  situation  of  four  identical  diodes 
and  where  other  circuit  parameters  are  balanced  regardless  of  the  actual 
parameter  values.  The  model  predictions  for  a set  of  four  identical 
diodes  is  shown  in  TABLE  5.  The  sets  of  measurements  in  TABLES  1 and  2 
do  not  approach  this  ideal  mixer  response  behavior. 

Another  characteristic  of  the  perfectly  balanced  response  calcula- 
tions is  that  the  remaining  (where  both  P and  Q are  odd,  TABLE  5)  res- 
ponses are  similar  in  level  to  their  values  for  the  slightly  unbalanced 
data  (e.g.,  TABLE  3).  Apparently,  these  responses  are  not  strongly 
dependent  on  the  degree  of  balance. 
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SECTION  5 

POWER  SERIES  COMPARED  WITH  TAYLOR- FOURIER  SERIES 

The  most  significant  differences  between  the  power  series  and 
Taylor-Fourier  series  with  regard  to  modeling  a receiver  mixer  are  as 
follows : 


1.  Using  the  power  series,  the  mixer  nonlinearity  is  repre- 
sented by  the  L0  amplitude  and  a series  of  coefficients,  the  minimum 
number  of  which  is  equal  to  the  highest  sum  of  the  harmonic  numbers  of 
the  L0  and  RF  signals  involved  in  the  mix  (i.e.,  P and  Q,  respectively). 

2.  Using  the  Taylor-Fourier  series,  representation  of  the 
mixer  nonlinearity  requires  a matrix  of  coefficients;  one  coefficient 
for  every  combination  of  P and  Q involved. 

Use  of  the  power  series  for  the  prediction  of  spurious  response 
or  intermodulation  provides  results  which  differ  substantially  from 
measured  results.  An  examination  of  two  of  these  erroneous  results 
follow. 


1.  For  a fixed  input  level,  VRp,  the  output  level  for  a P=i, 
Q=j  spurious  response  is  calculated  by  multiplying  the  P=j , Q=i  output 

level  by  the  factor 


<-1 * * * VL0/VRF') 


which  produces  results  contrary  to  measured  results. 

2.  For  an  increase  in  the  LO  voltage  by  a factor,  a given 
spurious  response  output  level  will  increase  by  that  factor  raised  to 
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the  p — power;  e.g.,  an  increase  of  the  LO  by  10  dB  will  increase  a 
P=5  response  by  SO  dB.  In  practice,  however,  a decrease  usually  occurs. 

Thus,  a receiver  mixer  model  based  on  the  power  series  will  generally 
lead  to  confusing  and  inaccurate  predictions  while  the  Taylor-Fourier 
series  models  are  consistent  with  experienced  results. 
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SECTION  6 
CONCLUSIONS 

1.  The  use  of  power  series  time-constant  coefficients  to  predict 
or  model  mixer  spurious  or  RIM  responses  should  be  avoided. 

2.  The  use  of  Taylor-Fourier  series  coefficients  to  predict  and 
model  mixer  spurious  or  RIM  responses  can  provide  reasonably  accurate 
results.  Models  are  described  for  two  types  of  mixers  and  instructions 
are  included  for  developing  additional  models. 

3.  Reasonable  estimates  of  high-order  RIM  responses  can  be  made 
on  the  basis  of  measured  values  of  low-order  RIM  responses.  This 
estimation  problem  is  more  difficult  for  situations  involving  double- 
balanced-diode  mixers  because  many  of  the  RIM  responses  depend  on  the 
degree  to  which  the  mixer  is  balanced.  A rule-of-thumb  factor  for 
dealing  with  this  situation  is  provided  in  APPENDIX  A. 
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APPENDIX  A 

COSAM  REQUIREMENTS  FOR  INCLUSION  OF  ADDITIONAL 
INTERMODULATION  (IM)  PRODUCTS 


BACKGROUND 


One  type  of  IM  problem  exists  for  which  measured  data  will  rarely 
be  available,  because  this  type  is  not  considered  in  MIL-STD-449  mea- 
surement procedures.  Many  cases  of  measured  IM  that  are  not  predicted 
in  a recent  evaluation  of  COSAM  performance  in  the  high  frequency  (HF) 
band9  could  be  explained  by  products  of  the  type  noted  in  Equation  A-l. 
COSAM  is  a computer  model  for  the  prediction  of  communications  systems 
performance  in  a cosite  configuration. 


The  Frequency  Relationships 

Recent  measurements  at  ECAC  that  involved  an  HF  receiver  have 
revealed  a type  of  IM  response;  this  response  has  been  ignored  in  the 
past  and  is  a subset  of  the  following  equation: 


:mfi  : nf2  : pfLo  ■ £if 

where 


(A-l) 


m , n,  p 


positive  integers 

input  signal  frequencies 

the  receiver  first  local  oscillator 
frequency 

the  receiver  first  IF. 


9M.  Lustgarten,  et  al,  Performance  of  the  Cosite  Analysis  Model  (COSAM) 
for  Selected  HF  Equipments , ESD-TR-76-010,  ECAC,  Annapolis,  MD,  Novem- 
ber 1976. 
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A 


Equation  A-l  represents  a set  of  equations. 

The  usual  IM  interaction  involves  the  following  subset  of  Equation 

A-l: 


+mfx  + nf2  = fQ  (A- 2) 

(i.e. , for  p = +1) 


where 


= the  receiver  tuned  frequency  = 


Note  that  if  n=0.  Equation  A-l  describes  a set  of  common  spurious 
responses. 

Another  type  of  interaction  noted  in  the  HF  measurements  was: 

+mf : + nf 2 = fJF  (A- 3) 

(i.e. , for  p = 0) 

Measurements  made  using  MIL-STD-449  procedures  only  consider  interac- 
tions that  involve  Equation  A-2  and,  more  recently,  a subset  of  Equa- 
tion A-3  (i.e.,  the  subset  where  m = n = 1). 

The  remaining  subsets  of  Equation  A-l  (i.e.,  for  p = +1  and  p > 2) 
were  also  noted  in  the  HF  measurements.  The  IM  interactions  normally 
considered  are  subsets  of  Equation  A-l  such  that  p = +1  (the  sign  is 
determined  by  receiver  design)  and  m + n = 2,  3,  5 and  7.  Recent  mea- 
surements10 indicate  that  these  subsets  do  not  take  into  account  approxi- 
mately 80%  of  the  potential  (2-signal)  IM  cases.  Most  of  these  additional 
IM  cases  were  for  p = 0,  +1  and  p = 2 and  3.  A few  were  for  p = +1  and 

10ECAC  lab  memo:  "KWM-2  IM  Tests  in  Support  of  Model  Development," 

April  1976. 
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m + n = 4 and  6.  Some  study  and  additional  measurements  are  needed  to 
determine  the  most  significant  cases. 

Power  Transfer  Relationships 


Given  that  the  frequency  relationships  exist  for  the  usual  2- signal 
IM  interaction,  the  level  of  RIM  interference  can  be  computed  using: 

P.  = m (P  -3  ) + n (P, -g.  ) - K (or  k'  ) dBm  (A-4) 
ino  *■  a a'  b bJ  m,n  m,n'  v 

where 


P. 

mo 


P 

a 


K 


m,n 


K 


m,n 


the  equivalent  on-tune  input  interference  power  that 
results  in  the  same  output  interference  as  the  com- 
bined IM  effect  of  P and  P, 

a b 

the  receiver  input  level  due  to  source  a 

the  receiver  input  level  due  to  source  b 

the  receiver  rejection  prior  to  the  nonlinearity  at 
frequency  of  source  a 

the  receiver  rejection  prior  to  the  nonlinearity  at 
frequency  of  source  b 

the  normalized  nonlinear  conversion  coefficient  for 
order  m + n of  the  mixer 

the  normalized  nonlinear  conversion  coefficient  for 
the  RF  amplifier. 


Because  the  usual  IM  involves  only  the  p = +1  case,  only  two  sub- 
scripts of  K and  K'  were  needed.  Methods  were  required  to  calculate  K 
values  for  the  additional  types  of  IM  that  involve  interactions  occurring 
primarily  in  the  mixer  circuits  of  receivers.  Measured  data  will  rarely 
be  available. 
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The  NCAP  Computer  Program 

The  Nonlinear  Circuit  Analysis  Program  (NCAP),  a computer  program 
obtained  from  RADC,  permits  the  analyst  to  estimate  a class  of  nonlinear 
effects  in  receivers  that  involve  small  signals.  However,  mixers  are 
generally  operated  as  "large  signal"  nonlinear  devices;  i.e.,  local 
oscillator  levels  that  are  necessary  for  efficient  mixer  operation  are 
large  enough  such  that  basic  assumptions  in  the  NCAP  models  do  not  gen- 
erally apply. 

EXTRAPOLATION  OF  K VALUES 
m,n 


Predictable  trends  in  response  levels  exist  for  most  mixer  types 
and  should  greatly  ease  the  implementation  of  additional  IM  products 
for  COSAM.  These  are  evident  by  the  Taylor-Fourier  series  models  and 
are  supported  by  measured  data. 

The  COSAM  parameters,  Km  n>  are  directly  related  to  Taylor-Fourier 
coefficients,  a..  n for  Q = m + n.  In  general,  however,  the  values  for 

> X 

ap  are  necessary  to  predict  the  additional  IM  products. 

* > X 

The  antilogs  of  COSAM  parameters,  Km  , are  inversely  proportional11 

to  the  Taylor-Fourier  series  coefficients,  a,  . As  mentioned  earlier, 

l,m+n 

the  only  subset  of  Equation  A-l  presently  considered  by  COSAM  is  the  P=1 

subset  (note  Q = m + n) . Therefore,  Km  may  be  considered  a measure  of 

a.  The  objective  is  to  find  the  relationship  between  a,  and 

l,m+n  J ? l,m+n 

the  more  general  an  . For  a given  m + n,  the  a_,  are  the  Fourier 
P,m+n  P,m+n 

series  coefficients  of  the  periodic  time  functions  am+nCt)* 

11Maiuzzo,  M.,  "Models  of  Nonlinear  Receiver  Interactions,"  1975  IEEE 
EMC  Symposium  Record,  San  Antonio,  TX,  October  1975. 
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To  extend  the  COSAM  capability,  K values  for  other  than  the  P=1 

m,  n 

case  are  needed.  For  purposes  of  discussion,  the  notation  K will 

def  m> n>  CP) 

be  employed  where  the  usual  K = K from  the  proportion  (see  Ref- 

r m,n  m,n(lj  r r 

erence  11): 


K = K + 20  log  (a.  /an  ) 

m,n(P)  m,n(l)  6 v l,m+n  P,  m+n' 

For  the  cases  of  most  interest,  i.e.  where  m + n > 2, 
is  usually  less  than  a few  dB. 


(A- 5) 

the  second  term 


Some  exceptions  may  be  noted  in  the  case  of  balanced  mixers.  For 
m + n odd  (the  usual  IM  case  considered  by  COSAM)  and  P odd,  the  second 
term  of  Equation  A- 5 can  be  shown  to  have  approximately  a value  of  -20 
log  P.  For  m + n odd  and  P even,  the  value  of  the  second  term  of  Equa- 
tion A-S  is  dependent  on  the  degree  of  balance  of  the  mixer  (-20  dB 
seems  typical).  The  other  case  considered  by  COSAM  is  m + n = 2,  where 
the  second  term  of  Equation  A-5  does  not  seem  to  have  a simple  relation- 
ship with  the  value  of  P. 

BASIS  FOR  EXTRAPOLATION 


The  ap  ^ are  the  Fourier  series  coefficients  of  the  periodic  Taylor 
series  coefficients  a^Ct).  Given  the  general  shape  of  an  a^t),  reason- 
able estimates  of  its  frequency  spectrum  or  Fourier  series  can  be  made. 

Most  mixers  behave  as  switches,  either  on-off  or  plus-minus  in  the 
case  of  balanced  mixers.  Thus  the  instantaneous  transconductance,  a^(t), 
can  be  approximated  by  a trapeziodal- like  waveform.  The  IM  products  of 
interest,  however,  involve  the  Q > 2 coefficients;  e.g.,  a2(t)  is  one-half 
the  rate  of  change  of  with  respect  to  the  RF  driving  voltage,  as  a func- 
tion of  time.  This  coefficient  has  its  maximum  value  when  the  mixer  is 
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turning  on  or  off  (a  much  shorter  period  of  time  than  the  width  of  a^t]). 
Thus,  the  spectrum  of  a^Ct)  approaches  the  spectrum  of  an  impulse,  and 
stays  fairly  constant  with  P . The  higher-Q-order  aq(t)  coefficients  are 
more  "impulse-like."  Note  that  they  only  appear  as  impulses  over  the 
range  of  P values  of  interest.  Given  the  actual  pulse  widths  of  the 
aq(t),  bandwidths  may  be  calculated.  In  the  case  of  a perfectly  balanced 
mixer,  the  even-order  aq(t)  coefficients  are  always  zero;  as  a result 
the  ap  q values  are  also  zero  for  even  values  of  Q.  However,  20  dB  addi- 
tional rejection  seems  to  be  typically  realized.  Also,  for  perfectly 
balanced  mixers,  with  Q odd  and  P even,  the  ap  n coefficients  are  non- 
existent  because  of  odd  symmetry  of  aq(t)  about  the  time  where  the  L0 
changes  sign.  For  these  coefficients,  20  dB  more  rejection  than  the 
(P-1,  Q)  response  is  typically  achieved 


COEFFICIENTS  INVOLVING  Q=1 


As  mentioned  above,  the  a^(t)  waveform  is  almost  trapeziodal.  Thus, 
the  coefficient  values  a p ^ have  an  amplitude  envelope  which  can  be  cal- 
culated using: 


Mp  = 1/P, 

npf  _ t <i 
LO  r 

Mp  = 1/P? 

npf,.t  >i 
LO  r- 

(A-6) 


where 


Mp  = calculated  envelope  of  ap  ^ coefficients.  The  input 
power  for  a standard  response  for  a (P,l)  response 
will  be  -20  log  Mp,  in  dB,  relative  to  the  (1,1)  input 
power 

tr  = rise  time  of  the  a^(t)  waveform,  approximated  by  a 
trapeziodal  waveform. 

These  near-impulses  are  nearly  90°  out  of  phase  with  the  LO;  therefore 
the  coefficients  will  increase  in  value  with,  and  in  some  cases  may 
even  be  proportional  to  P. 


I 
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This  envelope  is  plotted  in  Figure  A-l  for  one  of  the  measured  DBDM's. 

The  rise  time  was  taken  from  the  computer  program  output.  Note  that 
the  peaks  of  both  the  measured  and  the  predicted  values  are,  with  one 
exception,  near  the  envelope.  However,  several  predicted  values  lie 
in  nulls  that  are  deeper  than  those  for  the  measured  responses.  The 
cause  for  these  optimistic  DBDM  predictions  is  not  precisely  known. 

The  discrepancies  are  the  result  of  an  imbalance  in  the  circuits.  It 
is  more  realistic  to  employ,  for  ap  A predictions  involving  DBDM's,  the 
envelope  of  the  Q=1  data  and  an  additional  10  or  20  dB  rejection  for  even 
order  P values. 


Q = harmonic  number  of  input  signal  frequency 
P = harmonic  number  of  LO  frequency 
(see  Equation  C-2  of  APPENDIX  C) 

Figure  A-l.  Measured  and  predicted  spurious  response 
values  for  DBDM. 
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APPENDIX  B 

APPLYING  MIXER  SR  REJECTION  VALUES  TO 
RECEIVER  REJECTION  CALCULATIONS 


BACKGROUND 


When  estimating  receiver  rejection  to  a spurious  response  the  mixer 
SR  rejection  values  should  not  be  used  directly,  and  the  rejection  (dB) 
of  the  receiver  front-end  added  to  the  mixer  rejection  is  not  sufficient 
in  most  cases.  The  latter  method  has  been  used  in  the  past. 

The  mixer  rejection  to  a given  (P,Q)  spurious  response  R n is  the 

* > X 

difference  between  the  input  levels  (in  dBm)  for  the  (1,1)  response  and 
for  the  (P,Q)  response.  A standard  output  response  level  is  fixed  in 
both  cases;  usually  this  output  level  is  slightly  greater  (say  10  dB) 
than  the  noise  level  at  the  output.  However,  because  of  the  nonlinear- 
ity involved,  the  amount  of  rejection  is  highly  dependent  on  the  stan- 
dard response  level  chosen  for  the  test.  Moreover,  in  a receiver  the 
noise  levels  at  the  mixer  output  are  likely  to  be  quite  different  from 
trn  present  during  the  breadboard  test  of  an  isolated  mixer  circuit. 
Fort  <nately,  a simple  correction  is  possible  if  the  amount  (E)  is  known 
by  which  the  output  standard  response  level  differs  from  breadboard  to 
working  receiver. 

A SIMPLE  CORRECTION 


The  receiver's  rejection  can  be  calculated  using: 

eSR  = RP,Q  + eFE  - M-V  E/Q  (B-l) 
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3 


SR 


P.Q 

E 


“M 


FE 


the  receiver's  rejection  to  the  spurious  response  in 
question,  in  dB.  That  is,  the  difference  between  the 
receiver  input  level  required  to  cause  a receiver 
standard  response  of  (usually)  10  dB  (S+N)/N,  first  at 
the  receiver's  tuned  frequency  and,  second,  at  the 
spurious  response  frequency. 

the  mixer's  rejection  to  that  SR,  in  dB 


the  level  at  the  mixer  output  (if  the  receiver  input 
is  on-tune)  that  results  in  a receiver  standard  res- 
ponse (i.e.,  the  receiver  input  is  at  the  receiver 
sensitivity  level),  in  dBm 

the  mixer  output  reference  level  used  during  the  SR 
rejection  measurements  of  the  mixer  alone,  in  dBm 

the  rejection  relative  to  on- tune  provided  by  the 
receiver's  front  end  stages  to  a signal  at  the  SR 
frequency,  in  dB. 


AN  EXAMPLE 

In  TABLE  2 the  rejection  to  a (2,2-)  SR  is  82  dB  for  a HP  model 
10514  double-balanced  diode  mixer  (DBDM) . Here  R = 82  dB.  The  out- 
put  reference  level  is  -111  - 8 = -119  dBm;  the  conversion  loss  for  the 
(1,1+)  response  was  measured  at  8 dB.  Thus,  EM  = -119  dBm.  Now,  if 
this  mixer  is  used  in  a receiver  where  the  front-end  gain  is  24  dB  and 
whose  sensitivity  is  -105  dBm,  the  mixer  output  level  E will  be  -105  + 
24  - 8 = -89  dBm.  The  parameter  E in  Equation  B-l  then,  E = -89  - 
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(-119)  = 30  dB.  (Note  also  that  the  conversion  loss  of  8 dB  is  can- 
celled in  the  computation  of  E and  its  value  is  arbitrary  for  these  cal- 
culations.) Thus,  from  Equation  B-l,  the  receiver's  rejection  to  the 
(2,2,-)  response  is: 

BSR  = 82  + epE  - (2-1)  (30)/2 

or 

esR  = 67  + Bfe,  in  dB  (B-2) 

This  example  is  presented  graphically  in  Figure  B-l. 


Figure  B-l.  Graphical  presentation  of  example  problem. 
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AN  ADDITIONAL  CONSIDERATION 


Harmonic  levels  generated  by  the  radio  frequency  (RF)  amplifiers 
can  be  an  important  factor  in  a receiver's  rejection  to  a spurious 
response  (SR),  or  IM  response  for  that  matter.  For  example,  refer 
to  TABLE  2.  The  rejection  to  the  P=2,  Q=1  (+)  response  is  44  dB  less 
than  the  rejection  to  the  P=2,  Q=2  (+)  response.  Because  of  its  fre- 
quency, the  second  harmonic  of  the  signal  entering  the  receiver  at 
the  (2,2,+)  SR  frequency  will  pass  through  the  mixer  with  the  rejec- 
tion accorded  to  a (2,1,+)  response.  It  would  not  be  difficult  to 
construct  an  example  whereby  this  mechanism  would  determine  the 
receiver's  rejection  to  a SR.  It  would  be  important  to  include  this 
mechanism  if  the  SR  frequency  fell  below  the  tuned  frequency  and/or 
if  the  RF  amplifier  interstages  were  extremely  selective.  In  our 
example,  it  would  only  be  necessary  for  the  second  harmonic  of  the 
incoming  signal  to  be  44  dB  below  its  fundamental  (upon  reaching  the 
mixer)  in  order  for  that  mechanism  to  dominate  the  receiver's  response. 

Additional  stray  signals  may  also  be  introduced  into  the  mixer 
from  the  base  or  master  oscillator  of  from  other  circuits  in  the 
receiver,  although  generally  at  much  lower  levels. 
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APPENDIX  C 

"A  COMPUTER  PROGRAM  FOR  THE  PREDICTION  OF  DIODE  MIXER  SPURIOUS 
RESPONSE  LEVELS  FOR  CONDITIONS  INCLUDING  LARGE  LOCAL  OSCILLATOR  LEVELS" 

Mike  Maiuzzo 


Reprint  as  seen  in  IEEE  International  Symposium  on  Electromagnetic 
Compatibility , Seattle,  WA,  2-4  August  1977. 
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A COMPUTER  PROGRAM  FOR  TOE  PREDICTION  OF  DIODE  MIXER  SPURIOUS 
RESPONSE  LEVELS  FOR  CONDITIONS  INCLUDING  LARGE  LOCAL  OSCILLATOR  LEVELS 


M.  A.  Maiuzzo 
IIT  Research  Institute 
ECAC,  North  Severn 
Annapolis,  MD  21402 


Summary 


The  conversion  transconductance  of  mixers  in- 
creases with  increasing  local  oscillator  drive  level. 
The  rejection  of  unwanted  nonlinear  outputs  such  as 
intermodulation  and  spurious  responses  also  generally 
improves  with  local  oscillator  level.  Small-signal 
analysis  techniques,  including  conventional  employment 
of  the  Taylor  series  and  Volterra  series,  are  not 
applicable  under  these  conditions  and  can  lead  to 
large  errors  in  the  estimation  of  rejection  values. 

A procedure  for  the  calculation  of  wanted  and 
unwanted  response  levels  is  presented.  The  procedure 
employs  techniques  applicable  to  conditions  of  either 
large  or  small  local  oscillator  drive  levels.  Accurate 
predictions  are  shown  to  be  obtainable  based  on  know- 
ledge of  linear  circuit  parameters,  direct -current 
diode  characteristics  and  applied  voltages.  Compar- 
isons of  predictions  with  measured  nonlinear  response 
levels  are  made  and  the  method  is  shown  a be  accurate. 

The  model  also  indicates  that,  for  practical 
mixers,  the  responses  are  not  overly  sensitive  to  the 
diode  parameters.  Good  prediction  accuracies  are  ob- 
tained with  uncertainties  of  up  to  2S%  in  the  diode 
parameter  a and  a factor  6f  ten  in  the  reverse  leakage 
current. 


process  which  provides  the  Fourier  coefficients  also 
provides  the  smoothing  of  the  predictions  with  resp«;. 
to  the  input  data. 

Teman6  has  described  the  technique  in  a dis- 
cussion of  spurious  responses.  His  example  consists 
of  the  Fourier  expansion  of  the  transconductance,  the 
first  term  of  the  power  series.  RADC/SIGNATRON7  dis- 
cuss the  technique  in  more  detail.  A spurious  re- 
sponse is  said  to  occur  when  an  input  signal  is  at  a 
frequency,  fgp,  such  chat  equation  (2)  is  satisfied. 


PfL0  ifIF 


SP 


(2) 


where 

f^Q  » LO  frequency 

fjp  » the  IF,  or  desired  output  frequency 

P » the  harmonic  number  of  the  LO  frequency 
contributing  to  the  response  at  fjp 

0 * the  harmonic  number  of  the  input  signal 
frequency,  fgp,  contributing  to  the 

response  at  flp. 


Introduction 


Numerous  investigators  have  developed  models  of 
diode  mixer  nonlinear  responses.  However,  reported 
procedures  have  in  general  been  unsatisfactory  for 
practical  predictions  of  spurious  responses,  based 
soieiy  on  diode  DC  characteristics.  In  realistic  sit- 
uations, the  local  oscillator  (LO)  level  is  large.  The 
analyses  of  references  1 and  2,  for  example,  rely  on 
the  small-signal  assumption. 

One  technique  which  has  found  success  in  large- 
signal  intermodulation  analyses  is  the  Fourier-Bessel 
series  expansion. J It  has  recently  been  applied  to 
the  analysis  of  intermodulation  in  Classes  3 and  C 
transmitter  amplifiers . 4 »* 

In  this  paper,  a Taylor-Fourier  series  expansion 
is  applied  to  the  prediction  of  diode  mixer  spurious 
responses.  Tne  equation  for  the  static  diode  transfer 
is  given  byl : 

i - i0  - l]  (1) 


where 

i « diode  current 
i0  ■ reverse  leakage  current 

v • diode  voltage 
R^  • diode  bulk  resistance 

The  model  is  based  on  the  following  concept.  Each 
power  series  coefficient  of  the  expansion  of  the  diode 
current  is  treated  as  a periodic  time  variable.  The 
amplitude  and  period  of  each  is  determined  by  the 
amplitude  and  period  of  the  LO  drive.  Each  coefficient 
is  then  expanded  in  a Fourier  series.  The  integration 


It  is  shown  that  the  m ethod  of  analysis  reporter 
on, this  paper  does  allow  accurate  prcdictio.i  of  mixer 
spurious  responses. 


Summary  of  the  Approach 


The  static  DC  output  current  vs.  applied  source 
voltage  characteristics  is  expanded  in  a truncated 
Taylor  series: 


i ■ 


u-o 


(3) 


where 

• total  Thevenin  equivalent  applied  voltage 


II  * — 

\ q: 


SvQ 


Ci) 


and  is  evaluated  for  equal  to  the  "bias",  or 

operating  point  voltage.  For  this  case  the  circuit's 
DC  transfer  characteristic  is  given  by: 


*iR  j ) 


(3) 


where 

*T  * RS  * "b  * *1  ’ (6) 

Rg  • Thevenin  equivalent  source  resistance 

R^  » load  resistance 
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V 


(7) 


Noting: 

(vt-iR.r) 

c ‘ * i * 1 (Sa) 

i 

o 


Then: 

Li  = 1 

a J t 1 ♦ aUTCi  ♦ ioJ 

From  (4)  : 


aLi  ^ , 

1 *c*  R-|-(I  • iQ) 


C8) 


(9) 


■here  I is  the  value  of  i which  satisfies  equation  (S) 
with  the  bias  voltage  substitute  for  v^.  The  numerical 

solution  of  (5)  employs  both  Newton-Raphson  and  suc- 
cessive substitution  methods.  The  expressions  for  a, 
through  ij  .ire  presented  in  APPENDIX  I. 

Up  to  this  point,  the  analysis  has  been  the  same 
as  a standard  small  signal  analysis.  The  departure 
comes  from  treating  the  local  oscillator  drive  voltage 
as  a time-varying  bias  voltage. 


Figure  1.  IN82A  Diode  Current  and  Mixer  Trans- 
conductance  Computed  from  Equations 
(5)  and  (9). 


rewriting  equation  (3) ; 

7 

1 * £ a0(t)v  I (12) 

Because  aQ(t)  are  periodic  functions  of  tine,  each 
may  be  expanded  in  a truncated  Fourier  series. 

aP,y  t:OS  (13) 


Ur‘: 


ag(t) 


10 

E 

P=0 


10 

• E 

P=0 

7 

E 

Q-u 


As  a i esul  i , 


all  the 


r -r  — ‘ ■= 


time.  In  most  of  the  subsequent  discussion,  these  will 
be  referred  to  as  a^Ct).  Figure  1 illustrates  this 

by  an  example.  Both  I and  a^Ct)  have  the  same  period 


as  the  LG  waveform.  In  this  example,  moreover,  the  LO 
waveform  is  cosinusoidal;  the  bias  term  B in  equation 
(3)  is: 


B * VB  * VACC0S  9L0 


(10a) 


where 

Vg  is  the  DC  bias  voltage 

VAC  is  the  peak  voltage  of  the  Thevenin  equivalent 
LO  source  voltage 

9,0  has  vaiues  from  -n  to  -r  and  is  the  instantan- 
eous local  oscillator  phase  angle  (*7irf^jt). 

Thus,  B is  shown  to  be  periodic.  Note  also  that 


v 


t 


B » v 


RF 


01) 


where 

Vjy.  is  the  input  voltage  waveform  at  the  RF. 


First  w«  see  that  the  current  has  siauiy  fre- 
quency components.  In  general,  only  one  will  be 
of  concern,  the  one  appearing  at  the  IF  (i.e., 
f.p)  . For  example,  let  the  IF  = 3 MHz  and  the  LO 
frequency  (f,D)  be  33  MH2 . Then  any  signal  at  a 
frequency  such  that:  ’ 


{ 


RF 


35P*  3 


Q 


(IS) 


will  result  in  a component  at  3 MHz.  Moreover,  a 
signal  at  either  32  or  38  MHz  is  termed  a desired 
(or  image)  frequency;  these  correspond  to  P=Q=1. 
Signal  frequencies  of  67  and  73  MHz  will  also  re- 
sult in  3 MHz  responses  (P»2,  Q*l),  and  so  forth. 


Let  the  input  signal,  in  any  case,  be  given 

by: 


VKF  * VRF  COs2!IfUFt 


(16) 


The  peak  current  component  at  the  IF  is, 
14)  and  (16): 


1I’,Q  * ('$v RF*V ,Q 


. r>o 


from 

(17) 
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The  above  equation  must  now  be  related  to  the 
measured  quantities.  For  the  circuit  of  Figure  2, 
the  following  may  be  obtained. 

l’L  • 10  log  (R^  i2pQ)  dlia  CIS) 

Pav  * 10  106  W * 21  Mm  (19) 

«S 

where  P is  the  spurious  power  delivered  to  R 
and  PAV  is  the  power  available  from  the  source 

The  desired  response  is  defined  here  as  the  P«l, 
Q*1  response.  The  input  "on-tune"  power  level  POT 
required  to  yield  a given  IF  output  lev/sl  S (which 
is  defined  as  the  standard  response  output  level) 
is  given  by: 

P0T  * S ' 10  l0K  C»lKs)  * 20  log  aj.j  dBm  (20) 


The  input  available  power  level,  at  an  RF  fre- 
quency picked  to  result  in  a particular  spurious 
(P,3)  response,  called  P , may  be  related  to 

9 apux 

or 


p 

spur 


10  log  R_  ♦ 


/Q 


10  log  «s  ♦ 20  lug 

(r>o)  dBo  (21) 


sistance  is  printed,  along  with  the  second  harmonic 
current.  The  relative  levels  suggest  that  the  diode 
is  being  driven  extremely  hard.  This  is  the  approp- 
riate condition.  On  the  next  line  we  have  the  RF 
circuit  resistance.  From  equation  (14)  we  see  that 
this  is  1/a  , ohms,  the  diode  resistance  at  the  de- 

sired RF  muirthen  be: 

Rd  * C1/a0,l>  _ "t  0,ubs  l22) 

Finally,  we  have  the  matrix  of  input  (available) 
power  levels  required  to  produce  a standard  IF  re- 
sponse, in  this  case,  -116. 8 dBm. 
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Figure  3.  IMSPUR  Computer  Program  Sample 
Outpue. 

Comparisons  With  Measured  Data 


Figure  2.  Mixer  Circuit  Diagram 
Sample  Output 


Following  are  comparisons  with  measured  data. 

It  should  be  emphasized  that  the  predictions  were 
based  solely  on  the  DC  characteristics  of  the  diode 
and  a linear  circuit  analysis.  I.MSP'JR  is  not  an 
empirical  or  statistical  model  based  on  spurious  re- 
sponse or  dynamic  response  measured  data. 

The  measured  data  for  a diode  mixer  test  circuit 
are  shown  on  TABLE  1(a)  for  an  L0  voltage  (Thevenin 
equivalent)  of  0.468v,  peak.  The  corresponding 
IMSPUR  predictions  are  shown  in  TABLE  1(b).  Note 
that  forty-two  of  Che  predicted  values  (i.e.,  well 
over  half)  differ  by  less  than  1 dB  from  the  measured 
values.  The  average  difference  is  1 d3.  The  stan- 
dard deviation  of  the  differences  (or  error,  if  you 
like)  is  3.1  dB.  For  illustration,  some  of  these 
data  are  shown  on  Figure  4. 


A sample  output  sheet  of  the  IMSPUR  program 
is  presented  in  Figure  3.  The  inputs  are  printed 
first. 

On  the  next  line  it  may  be  seen  that  the  cir- 
cuit, although  extremely  hard-driven  is  not  very 
efficient.  (Note  that  the  desired  input  power  is 
user  supplied) . Design  features  to  achieve  more 
efficient  conversion  have  not  been  attempted.  Tech- 
niques could  include:  faster  L0  rise  time,  trap 
filters  at  the  LO  frequency  and  RF  to  by-pass  the 
load  resistor  (and  thus  increase  the  effective  diode 
voltages  at  these  frequencies)  and  matching  the 
sources  to  the  diode  resistance,  R.,  discussed  be- 
low. Next  are  some  calculated  parameters  of  inter- 
est. The  fundamental  frequency  local  oscillator 
current  estimated  to  be  present  in  the  load  re- 


A sensitivity  analysis  of  the  IMSPUR  prediction 
accuracy  was  performed,  as  a function  of  input  para- 
meter values.  The  results  are  presented  in  TABLES 
2(a)  and  2(b).  It  is  seen  that  the  standard  error  is 
minimum  when  the  values  of  i and  a arrived  at  through 
examination  of  the  DC  charac?eristics  are  employed. 

One  purpose  of  the  measurements  was  to  test  for 
a need  to  include  frequency  dependent  diode  parameters 
such  as  junction  capacitance.  For  this  purpose  refer 
to  TABLE  3.  In  this  table  are  the  input  signal  fre- 
quencies which  were  .applied  in  order  to  obtain  each 
(P,Q)  spurious  response.  There  are  two  values  for 
each  P,Q  combination;  note  the  t signs  of  equation 
(IS).  The  predictions  (TABLE  1)  associated  with 
signal  frequencies  above  300  MHz  are  noticeably  in- 
creasing in  error.  Although  accuracies  presently 
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Figure  4.  TABLE  1 Data  Plotted  for  Comparison, 
Through  Q-5 . 

obtainable  are  acceptable  for  many  purposes  in  the 
UHF  band,  there  seems  a need  for  additional  refine- 
ment here. 
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In  addition,  define  such  that 
v - RTaj  - 1 
then  we  obtain: 

2 


*2  * “3jU  /2 


(I-l) 

d-:j 

(i-jj 

d-4) 


a3  * aa2U  ^ * *J  r';3l 

a4  * * uK-r  (4ILjCiia,*»o;i  ,a.j 

♦ 8a,2  w)/12. 

aj  * ua^u^/5  ♦ alLj.  (30  ICja^a^a. 

* 24  * 72  a:a5  v*  2'ai  J.;  ^ ;> / 00 , (1-5) 

2 ■> 
a6  * °aS  v /b  * »rt(43-  H,.a2‘a. 
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V1!^  * 1B2  l;Ta,av‘.|  * 3I(m 


♦ 3840^11  * lCOa^ij.  i>)/3n0 

2 ■> 

a7  ' “o  **  /7  4 aKx  (3M0IUa,a_' 

♦ 28S0  Ryi^a^  ♦ 1440  ILj.a^a.a, 

♦ 1200  R^a,.-^  ♦ ESSOa.a^ 

♦ 2400a2a5  TCOa^y  )/252U 


(1-u) 


a--) 


The  measured  data  for  a condition  of  an  LO 
voltage  (Thevenm  equivalent)  of  1.6ov  peak  are 
shown  on  TABLE  4(a).  The  corresponding  IMSPUR  pre- 
diction* are  ™ TiBIF  4(b).  The  ts  aTe-com- 

parable  although  increased  error  is  seen  for  Q«5, 

6,  and  7 and  P*1  responses.  Two  possible  explana- 
tions are:  ..nability  to  handle  input  RF  levels 
approaching  the  LO  level  and,  second,  undetected  RF 
source  harmonic  leakage  into  the  test  circuit.  In 
any  case,  the  results  are  still  reasonably  accurate 
considering  the  one-hundred  dB  range  of  the  responses. 

CONCLUDING  COMMENTS 

While  additional  levels  of  complexity  could  have 
heen  employed  in  the  IMSPUR  computer  program,  the 
test  comparison  results  were  excellent.  One  can  ex- 
pect to  be  within  5 dB  of  the  mixer  sensitivity  to 
a given  spurious  frequency  for  17  orders  of  non- 
linearity. 

Additional  features  which  could  be  added  include 
consideration  of  the  diode  shunt  resistance  in  the 
Taylor  series  expansion. 

Reactive  elements  could  also  be  included;  note 
that  model  errors  increased  somewhat  for  frequencies 
above  300  MHz.  The  subject  is  covered  in  the  RADC 
report  in  more  detail. 

APPENDIX  I 

In  the  text,  the  periodic  time-varying  power- 
series  coefficient  a (t)  was  derived.  It  was  said 
that  a (t),  when  expanded  in  a Fourier  series, 
yielded  the  constant  coefficients  a ,,  for  p ■ 0,  1, 

2 Below  the  expressions  fer1  a (t)  for  Q • 

2,  3,  4,  . . , 7 are  presented.  For^convenience, 
these  will  revert  back  to  the  notation  used  previ- 
ously, i.e.,  a • a (t). 


Each  of  these  is  a periodic  tine  function  of  a 
periodic  time  function.  The  IMSPUR  computer  pro- 
gram expands  each  a,(t)  in  a Fourier  series  whose 
constant  coefficients  are  a„  n,  P«0,  1,  2,  . . . 

. . ..  10. 
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APPENDIX  D 

CALCULATING  DIODE  PARAMETERS 

TABLE  D-l  presents  dc  diode  measured  data  for  the  diode  used  in 
the  mixer  of  APPENDIX  C.  The  measurement  circuit  is  shown  in  Figure 
D-l.  For  each  dc  source  voltage  level,  the  voltages  at  points  (1) 
and  (2)  were  measured.  The  diode  voltage  is  the  difference  between 
these  two  data.  The  diode  current  is  the  voltage  at  point  (2)  divi- 
ded by  the  load  resistance  of  50  ohms. 

Examination  of  the  measured  data  clearly  indicates  a reverse  leak- 
age current  of  2 microamperes.  In  addition,  a diode  shunt  resistance 
of  125  kilohms  is  indicated;  note  that,  for  reverse  voltages  greater 
than  about  % volt,  the  reverse  voltage  has  a slope  of  about  125  kfi. 

We  now  have  two  of  the  four  diode  model  parameters.  By  submitting 
these  two  parameters  and  the  measured  data  for  forward  voltages  to  the 
computer  program  discussed  below,  the  other  two  diode  parameters 
are  obtained. 


Figure  D-l.  Measurement  circuit  for  dc  characteristic 
of  diode. 
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TABLE  D-l 

MEASURED  1N82A  DIODE  DC  CHARACTERISTIC 


REVERSE 

FORWARD 

V 

I 

V 

I 

(mV) 

(VA) 

(mV) 

(VA) 

-3988 

-263 

249 

2 

-2997 

- 67 

497 

6 

-1999 

- 23 

992 

16 

-1499 

- 12 

189 

231 

-1000 

- 6 

256 

926 

- 800 

- 6 

303 

2040 

- 600 

- 4 

341 

3340 

- 500 

- 2 

372 

4800 

- 400 

- 2 

404 

6320 

- 300 

- 2 

436 

7660 

- 200 

- 2 

465 

9150 

- 100 

- 2 

480 

10940 

- 50 

- 1 

620 

18520 

750 

26310 

970 

42730 

1180 

59360 
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Figure  D-2  provides  the  reader  with  a FORTRAN-V  listing  of  DCPAR, 
a computer  program  that  identifies  appropriate  values  of  the  diode  dc 
parameters  a (alpha)  and  R^  (bulk  resistance).  The  user  supplies  the 
following  data.  See  above  for  example  aquisition  of  the  first  two. 

1.  Diode  reverse  leakage  current,  i^ 

2.  Diode  model  shunt  (parallel)  resistance,  R^,  (results  are 
generally  very  insensitive  to  Rg^>100  kft.) 

3.  Pairs  of  diode  dc  current  and  forward  diode  vo.  ^es,  I 
and  V data,  respectively 

4.  The  number  (NVI)  of  current  and  voltage  data  pairs. 

A sample  output  is  provided  in  Figure  D-3.  The  first  column  is 
the  value  of  R^  selected.  The  second  column  is  the  range  of  calculated 
alpha  values  calculated  from  the  NVI  data  pairs  of  current-voltage.  The 
next  column  is  the  average  value  of  the  alphas  calculated  from  the  ND 
data  pairs  assuming  the  same  value  of  R,  . The  last  column  is  the  stan- 
dard deviation  of  alpha  values.  Since  alpha  and  R^  are  assumed  indepen- 
dent of  applied  voltage,  the  value  of  R^  that  yields  the  most  constant 
value  of  alpha  for  the  entire  set  of  submitted  data  pairs  is  seen  to  be 
R^  = 13  ohms  and  alpha  = 24.8  inverse  volts. 


I 
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C 

C FIND  ALPHA 

DIMENSION  V(100) » C ( lOO ) » ALPHA (100) 

READ  15*100)  NVI »RSH»CL 
100  FORMAT ( I2»2F8.6) 

DO  1 I=1»NVI 

READ  1 5 » 105 ) V < I ) ► C ( I ) 

105  FORMAT (2F8»6) 

1 continue 

RS=.l 

WRITE(6»110) 

110  FORMAT ( 16X t *RS*  »15X* ’MAX-MIN* , 15X» 
**AVE*  H5X»  * SIGMA*  » / » 16X > 2 ( 1H-)  » 
*15x»7(ih-) »15X»3(1H-) »15X»5(1h-) ) 

500  A=-1000. 

8=1000. 

TOTAL=0. 

DO  10  I=1»NVI 

ALPHA ( I ) =ALOG ( C ( I ) /CL -V < I) / ( CL*PSH) +1 ) 
*/ ( V ( I > * ( 1+RS/RSH ) -C ( I ) *RS ) 

A=AiVAX1  ( A#  ALPHA  ( I ) ) 

8=AMlNlO»ALPHA(I)  ) 

TOTAL=ALPHA ( I ) +TOTAL 
10  CONTINUE 
DIFF=A-S 

ave=total/nvi 

sigsq=g. 

DO  15  j=lrNVI 

15  SIG50=SIGSQ+ ( ALPHA < J) "AVE ) **2/NV I 
SIGMA=SQRT(SIGSQ) 

S»RITE(  6»  120 ) RS»OIFF»AVE»SIGMA 
120  FORMAT ( 10X » £1^*7 toX'El^.Ttex^El^^f 
*6X»E14.7) 

IF(RS-.l)  50 » 50 » 55 
50  R5=RS+.4 
GO  TO  500 

55  IF<RS-20.00QOGO>  oC»75»75 
60  RS=RS+.5 
GO  TO  500 
75  END 


Figure  D-2.  DCPAR  listing. 
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RS 

.1 000000+00 
.5000000+00 
.1000000+01 
.1500000+01 
•2000000+01 
.2500000+01 
.3000000+01 
.3500000+01 
.4000000+01 
.4500000+01 
.5000000+01 
.5500000+01 
.6000000+01 
.6500000+01 
.7000000+01 
.7500000+01 
.6000000+01 
.6500000+01 
.9000000+01 
.9500000+01 
•1000000+02 
.1050000+02 
•1100000+02 
•1150000+02 
•1200000+02 
•1250000+02 


•1400000+02 

•1450000+02 

.1500000+02 

.1550000+02 

•1600000+02 

.1650000+02 

•1700000+02 

•1750000+02 

.1800000+02 

•1850000+02 

.1900000+02 

.1950000+02 

•2000000+02 

•2050000+02 


aBest  a is  24.8 


■4 


MAX-MIN  A VE 


.1809652+02 

.1941942+02 

.1791671+02 

.1951967+0? 

.1768116+02 

.1964809+02 

.1743272+02 

.1978015+02 

.1717030+02 

.1991604+02 

.1689269+02 

.2005598+^2 

.1659855+02 

.2020021+02 

.1628636+02 

.2034900+02 

.1595443+02 

.2050264+02 

.1560081+02 

.2066144+02 

.1522333+02 

.2082576+02 

.1481949+02 

.2099598+02 

.1438646+02 

.2117254+02 

.1392095+02 

.2135592+02 

.1341919+02 

.2154667+02 

.1287679+02 

.2174539+02 

.1228860+02 

.2195279+02 

.1164861+02 

.2216966+02 

.1094966+02 

.2239691+^2 

.1018325+02 

.2263559+0? 

.9339121+01 

.2288694+02 

.8404858+01 

.2315239+02 

.7365237+01 

.2343364+02 

«62ol403+ol 

.2373275+02 

•5l69o52+ol 

.2405219+02 

.5168901+01 

.2439501+0? 

.5168750+01 

^24^6497+021 

.5448569+01 

.2516688+0? 

.7758108+01 

,2560693+n? 

.1049709+02 

.2609327+02 

.1379747+02 

.2663702+02 

.1785142+02 

.2725370+02 

.2295012+02 

.2796596+02 

.2955698+02 

.2680852+02 

.3845698+02 

.2983796+02 

.5109424+02 

.3115473+02 

.7044768+02 

.3295916+02 

.1038076+03 

.3572732+02 

.1749981+03 

.4097903+02 

.4326199+03 

.5805095+02 

.1555083+04 

-.6092563+02 

.3865443+03 

.1496420+02 

volts  *. 

Figure  D-3. 

DCPAR  output. 
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SIGMA 


.5254997+01 

.5192572+01 

.5111424+01 

.5026580+01 

.4937762+01 

.4644666+01 

.4746956+01 

.4644259+01 

.4536159+01 

.4422195+01 

.4301852+01 

.4174554+01 

.4039659+01 

.3896446+01 

.3744109+01 

.3581750+01 

.3408372+01 

.3222883+01 

.3024118+01 

.2810368+01 

.2582118+01 

.2337147+01 

.2076433+01 

•1803233+01 

.1527848+01 

•1278355+01 

a .1122547+g! 

•1172002+01 

.1489731+01 

.2043900+01 

.2803228+01 

.3778534+01 

.5021186+01 

.6630619+01 

.8784896+01 

.1181830+02 

.1642832+02 

.2433919+02 

.4124353+02 

.1029673+03 

.3638061+03 

.7871100+02 
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The  method  employed  in  the  program  makes  use  of  the  following  equa- 
tions . 

amperes  (D-l) 
volts  (D-2) 
amperes  CD- 3) 

Given  values  for  Rcu,  i~  and  R,  , then  a value  of  alpha  for  each 

on  U b 

pair  of  i and  v data  is  computed.  When  the  entire  set  of  NVI  pairs 
have  been  exhausted,  the  spread,  average>  and  standard  deviation  of 
resulting  alphas  is  printed.  The  value  of  R^  is  then  changed  and  the 
process  is  repeated.  The  increments  and  range  of  R^  are  fixed  in  the 
program.  DCPAR  input  formats  are  as  follows. 

CARD  1 


COLUMNS 

FORMAT 

DESCRIPTION 

1-2 

INTa 

NUMBER  OF  CURRENT  AND  VOLTAGE  DATA 

PAIRS,  NVI 

3-11 

FL.PT^ 

DIODE  MODEL  SHUNT  RESISTANCE,  R_„ 

on 

(OHMS) 

12-20 

FL.PT. 

DIODE  REVERSE  LEAKAGE  CURRENT,  l 

(AMPS) 

CURRENT  AND  VOLTAGE  DATA  CARDS 


COLUMNS 

FORMAT 

DESCRIPTION 

1-8 

FL.PT. 

DIODE 

VOLTAGE 

(VOLTS) 

9-16 

FL.PT. 

DIODE 

CURRENT 

(.AMPS) 

Integer. 

^FL.PT. - Floating  point. 
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APPENDIX  E 

NUMERICAL  SOLUTION  FOR  DOUBLE-BALANCED-DIODE  MIXER  CURRENTS 


The  circuit  diagram  of  the  double-balanced-diode  mixer  (DBDM)  is 
shown  in  Figure  6,  and  consists  of  four  independent  current  loops  or 
meshes.  In  the  computer  program,  the  applied  local  oscillator  is  a 
known  quantity  and  the  current  through  each  of  the  four  diodes  must  be 
calculated  from  Equation  16. 


If  the  four  diodes  were  each  replaced  with  a linear  resistor,  the 
problem  would  become  a trivial  linear  circuit  problem  that  could  be 
easily  solved  in  closed  form  with  simple  expressions  for  each  of  the 
unknown  currents;  these  currents  could  be  evaluated  for  each  possible 
applied  local  oscillator  voltage.  However,  the  presence  of  the  highly 
nonlinear  diodes  demands  an  alternative  approach. 


The  procedure  employed  is  a "relaxation"  method  in  which  initial 
assumptions  are  made  for  each  of  the  four  diode  currents.  The  value  of 
one  of  the  currents  is  then  "relaxed"  and  the  current  recomputed  with 
the  other  three  currents  regarded  as  fixed  at  their  presently  assumed 
values.  In  fact,  each  of  the  other  three  diodes  is  regarded  as  a cur- 
rent generator  supplying  the  assumed  current.  Under  these  conditions, 
the  equivalent  circuit  may  be  reduced  to  a single  mesh  containing  a 
diode,  an  equivalent  resistance,  and  an  equivalent  voltage  source.  The 
problem  then  becomes  one  of  solving  the  transcendental  Equation  E-l  for 

the  unknown  current  i . '. 

J 


i . 
3 


(v.-R.i.) 
3 3 y 


(E-l) 
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In  Equation  E-l,  the  terms  in  . and  aj  are  the  Ebers-Moll  parameters 


of  the  J 


th 


0,j 

diode;  is  the  total  resistance  of  the  single  mesh  equiva- 


lent circuit  corresponding  to  diode  j;  and  \T  is  the  equivalent  circuit 
voltage. 


Equation  E-l  is  the  Ebers-Moll  diode  model  (the  shunt  resistor  R^ 
is  ignored  here)  and  may  be  solved  for  the  unknown  current  by  means  of 
an  iterative  procedure.  To  insure  rapid  convergence  of  the  procedure 
by  which  the  diode  current  is  determined,  one  of  two  alternative  formu- 
lations is  selected  based  on  the  sign  of  the  voltage  . In  particular, 
for  V.  > 0,  the  Newton- Raphson  method  is  used  to  calculate  the  root  of 

gi<v  ■ ‘o.j  H [vvj  - Rjij)]  ■ *}  - v1  tE_2) 

while  for  < 0 the  root  of  G2 (i j ^ calculated  using: 


G,(i.)  = a . (v . - R.i.)  - In  (i./in  . + 1)  = 0 (E-3) 

2^  r Jv  3 J r 1 0,j  ; 1 

After  each  of  the  diode  currents  is  recomputed  in  this  fashion  it 
is  regarded  as  a known  quantity  while  each  of  the  other  three  currents 
is  "relaxed"  in  turn  and  subsequently  recomputed.  The  entire  process 
of  computing  all  four  currents  is  repeated  until  changes  in  successive 
estimates  of  each  of  the  currents  are  all  less  than  the  precision 

_9 

required  of  the  computation,  i.e.  10  amperes  in  the  actual  program. 

Experience  with  this  procedure  showed  that  in  general  the  process  con- 
-9 

verged  within  10  amperes  within  three  or  fewer  iterations  over  the 
set  of  four  diodes. 
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APPENDIX  F 
THE  CQ  EXPRESSIONS 


63 


ESD-TR-78-100 


Appendix  F 


ESD-TR- 78-100 


Appendix  F 


120  !ii  !fi 

J2  3e  3e2  3e3 

r s s s 


3i  V 32i 
r \ r 


s3irV 


J \3e3 
r s 


30  / 32lr\3  120  /9lr 

J2  W ) J 3 \3ec 


aW32ir 


, 3e  / \ 3e2 
s s 


r = 1,2, 3,4 


, 3i  36i 
7 r r 

J 3e  3e6 
r s s 

32i  3si 
21  r r 


' 3i  \ 2 35i 

r \ r 


' 3i  \3  34i 
r\  r 


J 3e2  3e5 

r s s 


-c„  3i  32i  34i 

250  r r r 


„ 33i  34i 

35  r r 


J2  3e  3e2  3e4 

r s s s 


3e  / 3e3 
'S'  s 


J 3e3  3e4 

r s s 


. / 3i  \2  32i  3 3i 

1,260  ( r \ r r 


3e  / 3e2  3e3 

1 s'  s s 


, 210  /fVf  f^fr  , 140  fir  /ffirV 
Jr  V3el  ) 3es  3es  '3eJ  / 


. 2520  / 3lr  t 32lr  + 2,520  /firV/fffr 

4 \3es/  3e,  Jj  , 

_ 630  fir  /fifr  f 720  /3lrV  r - 1,2, 3, 4 

Jf  3es  \,3ef  ) j6  \3eJ 


Note  that  the  units  of  are  that  of  inverse  volts  times  current  raised 
to  the  Qth  power  (i^/v) . 
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APPENDIX  G 

POWER  SERIES  DIVERGENCE  FOR  VACUUM  TUBES 

It  is  shown  that  the  power  series  for  vacuum  tube  amplifiers  and 
mixers  diverge  in  situations  of  interest.  A general  representation  of 
a tube  characteristic  (Chapter  6 of  Reference  3)  is  expanded  in  a power 
series.  The  ratio  test  for  absolute  convergence  is  applied,  and  the 
resulting  expression  is  discussed. 


The  tube  characteristic  is  calculated  using: 


3/2 

i = e 

P 

(G-l) 

where 

ip  = plate  current 

e = an  equivalent  plate  to  cathode  voltage. 

including  the 

effects  of  the  perveance,  grid  voltages, 
cation  factor,  as  applicable. 

and  amplifi- 

From  Equation  2,  a general  expression  for  the  n1^  order  power  series 
coefficient  is  given  by  Equation  G-2: 


3(-l)n(2n-5)l  Ip 
n!  (n-3) J 22n_3En 


n>3 


(G-2) 


where 


I 

P 

E 


quiescent  or  dc  plate  current 

quiescent  or  dc  equivalent  plate  to  cathode  voltage. 
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The  ratio  test  is  applied  by  finding  the  limit  as  n -*•  ® of  the 

absolute  value  of  the  ratio  of  the  (n+l)th  term  to  the  preceding  term. 

1 2 

If  the  limit  exceeds  unity,  the  series  diverges.  Define  this  limit 
as  r^,  then  using  this  and  Equation  G-2: 

r„  = lim  an+lCvRF5  =■  VRF  (G-3) 

Z , ,n  E 

n--  an  (vRp) 

where 

vRp  = the  ac  signal  voltage. 

As  an  example,  the  expression  for  E for  a triode  is  given  using: 

E = E ♦ Eb  (G-4) 

C - 

where 

Ec  = control  grid  to  cathode  bias  voltage 
E^  = plate  to  cathode  voltage 
U = plate  amplification  factor. 

Values  such  as  E^  = 100  and  u = 40  are  common  even  in  amplifiers. 
Typically,  E£  may  range  from  -1  to  -4  volts.  Therefore,  the  value  of 
E can  be  quite  small.  Consequently,  the  series  can  have  a very  small 
region  of  convergence. 


1 2 

Sherwood  and  Taylor,  Calculus,  Prentice-Hall  Inc.,  3rd  Ed.,  New  York, 
NY,  April  1954. 
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